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Abstract 
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to  estimate  solutions  of  equations,  in  particular  differential 
equations.  These  methods  can  be  used  to  obtain  information  on  the 
qualitative  behavior  of  solutions,  and  they  can  also  be  applied 
for  numerical  estimates.  This  report  provides  a  survey  on  the 
results  obtained. 


D '  s  t  r  V/  -  — - 1 

Avail  and/or  I 

!  Special 


A 


3 


Table  of  contents 

1  General  survey  4 

I.  OPERATOR  INEQUALITIES  7 

2  Inverse-positive  linear  operators  7 

3  Two-sided  bounds  for  second  order 

differential  operators  8 

4  An  estimation  theory,  range-domain 

implications  9 

5  Estimation  and  existence  theory  for 

vector-valued  differential  operators  11 

II.  SHAPE- INVARIANT  BOUNDS  AND  GENERALIZATIONS  13 

6  Description  of  the  operators  and 

estimates  considered  13 

7  Elliptic  operators  17 

8  Parabolic  operators  20 

III.  NUMERICAL  APPLICATIONS  26 

9.  Approximation  methods  and  a  posteriori 

error  estimates  for  two-point  boundary 

value  problems  26 

Literature  cited  30 

Appendix:  Contents  of  "Operator  Inequalities"  33 


4 


1.  General  survey 


The  research  reported  herein  is  concerned  with  estimates  for 
solutions  of  equations,  and  related  topics.  The  estimates  are  often 
combined  with  existence  statements.  A  series  of  applications  is  in¬ 
vestigated.  Particular  attention  is  given  to  numerical  applications. 

The  equations  considered  are  described  by  operators  M  of  the 
following  type: 

(1)  matrices,  and  nonlinear  functions  in  ]Rn  . 

(2)  Ordinary  differential  operators  of  the  second  order,  together 
with  boundary  operators. 

(3)  Ordinary  differential  operators  of  higher  order,  together 
with  boundary  operators. 

(4)  Ordinary  differential  operators  of  the  first  order,  together 
with  an  initial  condition. 

(5)  Vector-valued  ordinary  differential  operators  of  the  first 
and  second  order. 

(6)  Vector-valued  elliptic  differential  operators  of  the  second 
order,  together  with  boundary  operators. 

(7)  Vector-valued  parabolic  differential  operators,  together 
with  boundary  operators  and  initial  conditions. 

(8)  Abstract  differential  operators  of  the  first  order. 

(9)  More  general  abstract  operators. 


Thus  the  quantities  v  to  be  estimated  may  be  vectors,  real¬ 
valued  functions,  vector-valued  functions,  functions  with  values  in  an 
abstract  space,  or  more  generally,  elements  of  an  abstract  space.  The 
estimates  (inclusion  properties)  which  are  proved  state  that  v  be¬ 
longs  to  a  certain  set  K  .  In  particular,  we  consider  inclusion 
properties  such  that  v  £  K  is  equivalent  to  one  of  the  following 
relations  (a)  through  (e1). 


Estimates  for  abstract  elements  v  (with  applications  to  matric¬ 
es,  functions,  etc.) 


(a) 

< 

V  s  4. 

(two-sided  bounds)  . 

(a1  ) 

0 

V 

(positivity)  . 

(b) 

V 

e 

K 

(K  subset  of  an  abstract  space  with 
certain  properties)  . 
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(E) 


Estimates  for  vector-valued  functions  v 


(c) 

(d) 
(d1  ) 
(d" ) 


||  v  (x) ||  <  4.  (x) 
v(x)  €  <J,(x)G 
v  (x)  e  g 


(pointwise  norm  "bounds/  ||  ||  =  ^  . 

(shape-invariant  bounds,  4  (x)  6  3R1  ,  G  €  IRn)  j 
(invariance  statements)  . 


v(x)  £  4>^(x)G^  (£=1 , 2  ,  .  .  .  ,N)  (generalization  of  (d) )  . 


(In  these  estimates  the  variable  x  6  JRm  is  to  be  replaced  by 


(x,t)  €  3Rm+1  for  parabolic  problems.) 


Estimates  for  functions  v  with  values  in  a  Banach  space: 

(e)  v(t)  £  <J)(t)G  (G  a  subset  of  a  Banach  space  X  ,  in 

particular,  X  =  C§(fi))  . 

(e')  /  vT(x,t)  H  (x)  v(x,t)  dx  <  4. 2  ( t )  (special  case  of  (e)  ) 

Q 


The  basic  statements  proved  are  of  two  different  types  (U)  and 


(U)  Range-domain  implications: 

Mv  F  C  »♦  V  €  K  .  (1.1) 

Here  the  estimate  v  €  K  of  the  unknown  quantity  v  is  derived  from 
certain  properties  of  the  known  image  Mv  of  v  under  an  operator 
M  .  For  example,  if  M  is  a  differential  operator,  the  term  Mv  €  C  , 
in  general,  represents  a  differential  inequality,  or  a  set  of  differen¬ 
tial  inequalities.  Hence,  we  are  concerned  here,  in  particular,  with 
a  theory  of  differential  inequalities. 

Statements  of  the  form  (1.1)  can  be  used  to  obtain  estimates 
for  all  solutions  of  a  given  equation  Mv  =  r  : 

r  €  C  =*•  v  €  K  for  each  solution  of  Mv  =  r  (1.2) 

From  this  result,  one  can  often  derive  the  uniqueness  of  the 
solution.  This,  however,  is  not  the  main  object  in  proving  statements 
of  type  (U)  . 


(E)  Existence  and  inclusion  statements: 


r  £  C  -►  there  exists  a  solution  v  £  K  of  Mv  =  r  .  (1.3) 


In  addition,  modifications  of  (1.1)  are  considered,  where  certain 


properties  of  v  (such  as  positivity)  are  assumed  to  be  known.  Then 
we  have  an  implication  of  the  form 

(Mv  €  C  ,  v  €  K°)  -»  v  e  K  .  (1.4) 


Most  of  the  results  are  contained  in  the  monograph  "Operator 
Inequalities"  [1]  and  the  papers  [2, 3, 4,5].  Certain  numerical  re¬ 
sults  have  not  yet  been  published. 

This  report  is  divided  into  three  chapters.  Chapter  I  gives  a 
brief  survey  of  the  contents  of  "Operator  Inequalities".  This  book 
is  mainly  concerned  with  estimates  of  type  (a),  (a*),  (c) .  In  addi¬ 
tion,  methods  for  obtaining  results  on  more  general  estimates  (b) 
are  described.  The  operators  treated  are  abstract  operators  (9) , 
matrices  (1)  and  ordinary  differential  operators  (2)  through  (5). 
Chapter  II  of  this  report  describes  the  results  on  estimates  (d) , 
(d'),  (d"),  (e) ,  (e')  for  vector-valued  differential  operators  (5), 

(6) ,  (7)  and  abstract  differential  operators  (8) ,  obtained  in 
[2, 3,4,5].  Chapter  III  is  concerned  with  numerical  applications.  The 
appendix  yields  a  detailed  table  of  contents  of  [1], 
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I.  OPERATOR  INEQUALITIES  ! 

\ 

The  book  [1]  is  concerned  with  inequalities  that  are  described 
by  operators  or,  briefly,  with  operator  inequalities.  As  the  title 
suggests,  abstract  terms  are  used  in  developing  the  theory  and 
methods.  Abstract  results,  however,  are  not  considered  as  ends  in 
themselves,  but  as  means  to  obtain  results  for  concrete  problems. 

The  book  concentrates  on  matrices  and  (scalar-valued  and  vector¬ 
valued)  ordinary  differential  operators.  Because  of  size-limitations 
on  the  book,  partial  differential  operators  essentially  had  to  be 
omitted.  We  point  out,  however,  that  most  of  the  results  of  type 
(U)  for  ordinary  differential  operators  can  be  carried  over  without 
any  essential  difficulty  to  elliptic-parabolic  partial  differential 
operators  of  the  second  order. 


More  generally,  a  main  purpose  of  the  book  is  to  provide 
methods  of  various  kinds  which  can  also  be  used  for  problems  other 
than  those  treated  in  the  book. 


The  Sections  2  through  5  below  are  concerned  with  Chapters  II 
through  V  of  [1],  respectively. 

The  book  contains  more  than  330  references,  which  will  not  be 
listed  here. 


2.  Inverse-positive  linear  operators 


The  second  chapter  of  [1]  provides  a  unified  theory  of  linear 
operators  M  in  an  ordered  linear  space  R  such  that  for  v  €  R 

Mv  >  0  »•  v  >  0  (inverse-positivity) 

This  implication  is  of  type  (1.1)  where  Mv  €  C  is  equivalent  to 
Mv  >  0  .  The  basic  result  is  the  monotonicity  theorem  (Theorem  1.2). 

From  this  theory  most  of  the  known  results  on  M-matrices  and 
inverse-positive  second  order  differential  operators  are  derived,  as 
well  as  many  new  ones.  The  abstract  formulation  used  here  allows  one, 
for  example,  to  recognize  the  common  properties  of  these  two  differ¬ 
ent  operator  classes.  This  leads  to  a  theory  of  abstract  M-oper- 
ators.  Certain  general  methods  for  proving  inverse-positivity  are 
applied,  for  example,  to  differential  operators  of  higher  order. 

Inverse-positivity  (strict  inverse-positivity  etc.)  is  related 
to  a  series  of  other  theories,  and  it  can  be  used  in  many  applications. 
Here,  we  can  only  indicate  briefly  some  of  the  topics  which  are 
treated  in  [1]: 

Eigenvalue  theory  for  M-matrices,  second  order  differential 
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operators  and  M-operators,  in  particular,  the  Perron-Frobenius  theory 
and  generalizations. 

Convergence  theory  for  iterative  procedures  involving  matrices 
or  abstract  operators.  ~  ~  ~ 

Oscillation  theory  for  second  order  differential  operators. 

Boundary  maximum  principles. 

Convergence  theory  for  difference  methods  for  boundary  value 
problems. 

Error  estimation  for  approximate  solutions  of  initial  value 
problems  and  boundary  value  problems. 


Let  us  mention  some  typical  statements  of  the  corresponding 
theories,  which,  of  course,  hold  only  under  certain  assumptions 
which  we  cannot  formulate  here: 

"M  is  inverse-positive  if  and  only  if  each  real  eigenvalue  of 
M  is  positive"  (see  Corollary  2.14b  and  Theorems  3.24,  4.19). 

"The  iterative  procedure  for  Mv  =  r  converges  if  and  only 
if  M  is  inverse-positive"  (see  Theorems  2.19,  4.20). 

"M  on  [0,1]  is  inverse-positive  if  and  only  if  [0,1]  is  an 
interval  of  non-oscillation"  (see  Theorem  3.22). 

"The  difference  method  converges  if  the  differential  operator 
is  inverse-positive"  (see  Proposition  5.5). 

These  statements  show  that  each  of  the  sufficient  properties 
for  inverse-positivity  which  are  proved  in  Chapter  II  can  also  be 
used  in  other  theories  and  applications. 


3.  Two-sided  bounds  for  second  order  differential  operators 


Chapter  III  of  [  1  ]  treats  the  theory  of  two-sided  estimates 
<p  <  v  <  (j)  for  solutions  v  of  two-point  boundary  value  problems 
of  the  second  order.  (Initial  value  problems  of  the  first  order  are 
included  as  a  special  case.) 

First  we  describe  a  theory  on  inverse-monotone  differential 
operators  M  ,  which  have  the  property  that 

M<p  S  Mv  S  Mt|>  «*  ip  £  v  <  <J> 


Here,  <  denotes  a  po: itwise  inequalitiy,  and  M  is  an  operator  of 
the  form 
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-  a(x)u"(x)  +  f (x,u (x) ,u' (x) )  for  0  <  x  <  1 

Mu (x)  =-|  go(u{0),u'(0))  for  x  =  0 

gi  (u ( 1 ) , u*  ( 1 ) )  for  x  =  1 

Secondly,  we  treat  a  generalized  inverse-monotonicity ,  where 
the  pointwise  differential  inequalities  M<p  <  Mv  <  are  replaced 

by  weak  differential  inequalities.  It  is  shown,  for  example,  that 
differential  operator  which  is  monotone-definite,  i.e.  monotone  in 
the  sense  of  Browder-Minty,  is  also  inverse-monotone  (see  Theorem 
3.9). 

Furthermore,  Chapter  III  of  [  1  ]  deals  with  the  theory  of  comp 
ison  functions,  where  statements  of  the  following  form  are  proved: 


M<p  <  0  <  M4> 


there  exists  a  solution 
Mv  =  0  with  <p  <  v  <  i)> 


of 


and,  in  addition,  estimates  $  <  v'  <  V  for  the  derivative  of  the 
solution  are  obtained.  The  theory  is  derived  in  a  way  which  makes 
consequent  use  of  the  theory  of  inverse-monotone  operators. 

Finally,  the  theory  of  generalized  inverse-monotonicity  is  used| 
to  develop  a  theory  of  weak  comparison  functions ?where  the  pointwise 
differential  inequalities  Mq>  <  0  <  Mt|,  are  replaced  by  weak  differ 
ential  inequalities. 

The  examples  treated  include  singular  perturbation  problems, 
bifurcation  problems,  existence  proofs  for  positive  solutions  and 
numerical  error  estimates. 


In  developing  the  theory,  we  explain  (several)  methods  of  proof 
and  show  also,  how  these  methods  can  be  applied  to  other  problems, 
such  as  problems  with  periodic  boundary  conditions  and  certain  sin¬ 
gular  boundary  value  problems. 


4.  An  estimation  theory,  range-domain  implications 


Chapter  IV  of  l  1  3  describes  some  general  principles  which  can 
be  used  to  obtain  sufficient  conditions  for  statements  of  the  form 
(1.1): 

Mv  e  C  «*  v  £  K  ,  (4.1) 

where  C  and  K  are  subsets  of  certain  linear  spaces  R  and  S  , 
respectively.  The  results  can  be  used  to  obtain  sufficient  conditions 
on  M  such  that  (4.1)  holds  for  given  sets  C  and  K  .  They  may 
also  be  used,  however,  to  find  a  set  C  such  that  (4.1)  holds,  if 
M  and  K  are  given. 
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The  main  tool  is  a  continuity  principle  which  can  be  applied 
in  the  following  way.  One  constructs  a  family  of  sets  K 

A 

(O  <  A  <  y  <  co)  with  K0  =  K  and  *  K  for  A  >  0  such  that 

there  exists  a  minimal  A  with  v  £  K  and  that  for  A  >  0  the 

A 

element  v  belongs  to  the  boundary  (or  some  other  distinguished 
subset)  of  .  Then,  from  some  general  properties  of  the 

elements  v  £  ,  one  deduces  certain  relations  for  Mv  .  If  these 

relations  contradict  Mv  £  C  ,  the  statement  (4.1)  is  proved. 


This  continuity  principle  is  here  applied  to  obtain  results  on 
range-domain  implications  of  the  following  type: 

<  w  ;  (4.2) 


q>  <  V  <  (J;  (4.3) 


operators  in  ordered 


<  «>) 


(0  <  A  <  «>)  , 

respectively,  where  z\  and  z*  denote  elements  in  R  .  The 
results  on  inverse-monotonicity  are  applied  to  nonlinear  operators 
in  !Rn  and  vector-valued  functional  differential  operators  of  the 
second  order. 

For  initial  value  problems  of  the  first  order  a  second  approach 
is  described,  where  the  sets  K  are  defined  differently. 

A 

The  object  of  this  chapter  is  not  only  to  derive  sufficient 
conditions  for  special  statements  of  the  form  (4.1),  but  also  to 
provide  methods  which  can  be  used  for  operators  and  problems  not 
treated  here.  The  derivation  of  results  for  infinite  systems  of 
differential  equations  and  abstract  differential  equations  in  this 
chapter  are  to  be  considered  as  further  examples  for  the  application 
of  these  methods.  (See  also  Section  5.) 

Chapter  IV  of  [1]  contains  also  a  special  theory  on  statements 
of  the  form  (4.1)  for  linear  (and  concave)  operators.  This  theory 
is  applied  to  obtain  statements  of  the  form 


Inverse-monotonicity:  Mv  <  Mw  «*  v 

Two-sided  bounds  (with  Mu  =  H(u,u)) 

H  ( <p  ,h  )  <  Mv  <  H  (4>  ,h) 

for  all  h  with  <p  <  h  <  ^ 

In  the  theory  on  these  two  properties  for 
spaces  we  use  sets 

K^={u£R:u<w+  z^}  (0  <  X 

and 

K^  =  {u€R:  <p-z^<u<  +  +  z^} 


J. 
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-  4>  <  Mv  <  S'  -*  <p  <  v  <  4, 

for  second  order  linear  differential  operators  which  are  not  inverse¬ 
positive. 


5.  Estimation  and  existence  theory  for  vector-valued  differential 
operators  ‘ 


Chapter  V  of  [  1  ]  is  mainly  concerned  with  vector-va lued 
differential  operators  M  of  the  form 

f  -  a(x)u"(x)  +  b(x)u'(x)  +  f (x,u (x) ,u ' (x) )  for  0  <  x  <  1 


MU  (x)  =  K 

-  a°u '  (0)  +  f° (u  (0) ) 

for 

1 

| 

a  1u '  (1)  +  f1 (u ( 1 ) ) 

for 

with  u(x)  =  (u  j_  ( x ))  €  3Rn  and  diagonal  matrices  a(x)  ,  b(x)  ,  a0  , 
a1  .  Part  of  the  theory  is  also  formulated  for  functional-differen¬ 
tial  operators,  where,  for  example,  f  (x,u(x) ,u'  (x) )  is  replaced 
by  g(x,u(x),u'(x),u,u’)  .  Again,  operators  of  the  first  order  can 
often  be  treated  as  a  special  case. 

Operators  as  described  above  are  inverse-monotone  only  under 
very  restrictive  assumptions.  In  general,  one  requires  the  following 
two  conditions  on  the  components  ^  of  f  : 

(i)  M  is  weakly  coupled,  i.e.,  fj[  does  not  depend  on  any 
u’ (x)  with  k  *  i  . 

(ii)  M  is  quasi-monotone ,  i.e.,  fp  is  an  antitone  function 
of  each  uj^(x)  with  k  *  i  (antitone  =  non-increasing). 

If  both  these  conditions  are  satisfied,  the  theory  of  inverse- 
monotonicity  and  the  theory  of  comparison  functions  for  scalar¬ 
valued  operators  can  be  generalized  to  vector-valued  operators, 
without  major  modifications.  The  question  arises  what  kind  of  results 
can  be  obtained  without  requiring  both  these  conditions. 

Chapter  V  of  [  1  ]  contains  results  of  type  (U)  and  (E)  with 
estimates  of  the  form 


(a)  <p  <  v  <  4  ,  i.e.  , 

tp^(x)  <  v^(x)  <  4-p(x)  for  all  indices  i  and 

all  x  €  [0,1 ) 
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(c)  I i  v  ( x )  1 1  <  4(x)  for  0  <  x  5  1  , 

4  real-valued,  |j  |j  ,  'y*  ,  ,  y  an  inner  product  in  ]Rn. 

First,  these  estimates  are  derived  from  differential  inequali¬ 
ties  involving  v  (results  of  type  (U) ) .  In  the  case  (a)  one 
obtains  2n  coupled  differential  inequalities,  which,  in  addition, 
contain  two  parameters  h,k  £  ]Rn  .  For  example,  for  a  fixed 
x  €  (0,1)  ,  the  differential  inequality  for  4, ,  assumes  the  form 

(Mv),(x)  <  -  an(x)(i,!|(x)  +  b11(x)41(x)  +  f,(x,h,k)  .  (5.1) 

This  inequalitiy  has  to  be  satisfied  for  all  h,k  £  IRn  such  that 
h  =  (hi)  ,  h  1  =  4. 1  (x)  ,  <Pi(x)  <  hi  <  4>i  (x)  , 

k  =  (ki)  ,  ki  =  4 \ (x) 

The  occurrence  of  k  leads  one  to  require  assumption  (i)  ,  in  gener¬ 
al.  If  (i)  is  satisfied,  the  parameter  k  can  be  completely  elim¬ 
inated.  Then,  the  corresponding  theory  can  be  derived  from  the 
abstract  theory  on  property  (4.3).  Assumption  (ii)  need  not  be 
required.  If  (ii)  is  also  satisfied,  then  h  in  (5.1)  may  be 
replaced  by  4(x)  ,  so  that  no  parameter  occurs. 

For  estimates  (c)  one  obtains  a  differential  inequality  for 
4  which  again  contains  two  parameters  q,n  £  1Rn  : 

-  4"  (x)  +  <(q,q)4(x)  +  ,  f  (x,4  (x)  n ,  4’(x)n  +  4(x)q)^>  >  <(n,Mv(x)) 

for  0  <  x  <  1  ,  <(n.  n>  =  1  ,  <n,q>  =  O  . 

Here  none  of  the  assumptions  (i),  (ii)  need  to  be  required. 

There  are  other  essential  differences  between  the  theories  of 
the  estimates  (a)  and  (c) .  For  example,  if  f(x,y)  =  Cy  with  a 
matrix  C  ,  certain  eigenvalues  of  matrices  related  to  C  play 
an  important  role;  in  case  (a)  the  smallest  eigenvalue  of  the 

matrix  B  =  (b.  .)  ,  b .  .  =  c, .  ,  b.  .  =  -  |c,.|  ;  in  case  (c)  with 
T  'LJ“  ^  J  I 

<y,y>  =  y  y  the  smallest  eigenvalue  of  ^(C+CT)  . 


The  methods  used  to  prove  these  results  on  range-domain 
implications  can  also  be  applied  to  estimates  different  from  (a) , 
(c) . 


In  addition  to  these  results,  Chapter  V  yields  a  theory  on  the 
existence  of  solutions  v  which  satisfy  estimates  of  type  (a)  or 
(c) .  Here,  differential  inequalities  similar  to  those  described 
above  occur. 


* 
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II.  SHAPE- INVARIANT  BOUNDS  AND  GENERALIZATIONS 


In  this  chapter  we  describe  results  on  estimates  for  differential 
operators  which  were  presented  in  [2-5].  These  estimates  are  more 
general  than  estimates  by  two-sided  bounds  and  pointwise  norm  bounds. 
On  the  other  hand,  the  sufficient  conditions  derived  are  still 
comparibly  simple.  In  order  to  obtain  such  estimates,  one  essentially 
has  to  construct  one  (or  several)  real-valued  functions  satisfying 
certain  differential  inequalities. 


The  papers  [2]  and  [3]  treat  ordinary  differential  operators 
and  elliptic-parabolic  partial  differential  operators,  respectively. 
[4]  describes  a  different  approach  to  parabolic  operators,  which 
yields  somewhat  stronger  (and  additional)  statements.  [5]  is  a 
brief  survey  with  additional  results. 

Here,  we  shall  describe  the  basic  ideas  and  results  of  the 
theory  by  considering  certain  elliptic  and  parabolic  operators 
of  a  special  form.  The  papers  mentioned  treat  more  general  differen¬ 
tial  operators.  Moreover,  we  explained  in  [5],  that  the  results  can 
easily  be  generalized  to  certain  functional-dif f erential  operators. 
Furthermore,  the  methods  used  in  this  theory  may  also  be  applied 
in  proving  more  general  inclusion  properties. 


6.  Description  of  the  operators  and  estimates  considered 

In  this  report  on  the  results  in  [2, 3, 4, 5]  we  restrict  ourselves 
to  considering  certain  special  classes  of  elliptic  and  parabolic 
differential  operators. 

The  elliptic  operators  discussed  in  Section  7  have  the  form 

Mu(x)  =  -  D(x)  Au  (x)  +  f (x,u(x) ,ux(x) )  for  x  €  fi  (6.1) 

Mu  (x)  =  -  a  (x)  -|^(x)  +  B  (x)  u(x)  for  x  £  9fl  (6.2) 

where  fi  is  a  bounded  domain  in  IR™  ,  x  =  (x^)  ,  u(x)  =  (u^(x))  €  IRn 
ux(x)  =  (3u^/9xj^(x)  €  3Rn»m  ,  9/9v  denotes  the  interior  normal 
derivative  (or  some  other  directional  derivative  into  fi  ) , 

D  (x)  =  (dA(x)  &i  j )  €  3Rn»n  ,  d±(x)  £  0  ,  f  =  (fi)  :  fi  *  3Rn  *  3Rn'l"  -»  Kn 

a(x)  =  (0i(x)  6i  j )  €  ]Rn'n  ,  ai(x)  >0  ,  B(x)  =  (bik  (x) )  €  Kn  'n  , 

and  ]Rn'm  denotes  the  set  of  real  n*m  matrices. 
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In  the  following  sections  we  shall  often  use  a  simpler  notation, 
where  the  independent  variable  is  omitted  in  some  places.  For 
example,  we  may  write  (6.1)  in  the  form 

Mu  =  -  D  Au  +  f(x,u,u  )  on  the  set  ft 

A 

We  point  this  out  because  such  relations  will  be  combined  with 
side  conditions t where  the  independent  variable  cannot  be  omitted. 

For  operators  (6.1),  (6.2)  we  consider  estimates  of  the  form 

v(x)  €  <)<  (x)  G  for  x  £  ft  (6.3) 

v(x)  €  (x)  (£=1 , 2, . . . ,N)  for  x  £  ft  .  (6.4) 

In  (6.3),  G  denotes  a  given  closed  set  in  lRn  ,  v  £  R  =  Sn 
is  the  unknown  function  to  be  estimated,  £  S  a  function  to 
be  constructed,  and  S  =  Ci  (ft)  fl  Ca  (ft)  (for  simplicity)  . 

In  order  to  obtain  practicable  conditions,  we  assume  that  the 
set  G  c  ]Rn  is  described  by  a  real-valued  function  W  on  ]Rn  , 
such  that 


y  €  G  •*  W  (y )  <  1  . 

This  function  has  to  satisfy  certain  smoothness  conditions  and 
other  assumptions,  which  imply  that  G  is  star-shaped  with  respect 
to  0  and  that  G  has  a  smooth  boundary  r  =  {y  :  W  (y)  =  1 }  .  It 
is  also  required  that  the  normal  vector  W' (y)  at  a  boundary  point 
y  of  G  must  not  be  orthogonal  to  y  ,  i.e., 

ui(y)  :  =  W'  (y)y  >  0  for  W(y)  =  1  .  (6.5) 

The  estimate  (6.1)  is  equivalent  to 

V(v(x))  <  <j>(x)  for  x  6  ft  , 

where  V  is  the  Minkowski  functional  of  G  ,  which,  however,  need 
not  be  known  (see  [2]). 

For  constant  <)>{x)  s  <J;0  one  obtains  results  on  invariant  sets. 
For  variable  <1<  (x)  ,  the  size  of  the  set  4>(x)G  depends  on  x  , 
but  its  "shape"  remains  invariant.  We  speak  of  shape-invariant  bounds 

For  estimates  (6.3)  analogous  assumptions  are  made:  the  sets 
Gj,  are  described  by  functions  W£  ,  etc. 

The  parabolic  operators  M  considered  in  Section  8.1  have  the 


form 


Mu (x,t) 


(6.6) 


Mu (x, t ) 

Mu(x,0) 


ut(x,t)  -  D (x) 
-  a (x)  |^(x,t) 
u(x,0)  -  r (x) 


&u(x,t)  +  f (x,t,u(x) ,ux(x) ) 

for  x  £  ft  ,  0  <  t  <  T 
+  B (x)  u(x,t) 

for  x  £  ft  ,  0  <  t  <  T 
for  x  £  ft 


(6.7) 

(6.8) 


Here,  the  notation  is  essentially  the  same  as  above.  (In 

A  and  9/3v  are  differential  oper- 

form 

for  x£ft,0<t<T 

(£=1 , 2, . . . ,n)  for  x  €  ft  , 

0  <  t  <  T 

Here,  G  and  G£  are  sets  as  above,  v  £  R  =  Sn  ,  e  S  , 

4>£  €  S  ,  where  now  S  is  the  set  of  all  real-valued  functions  on 

n  x  [0,T]  which  are  continuously  differentiable  with  respect  to  t 
and  twice  continuously  differentiable  with  respect  to  x  on 

x  (0,T ]  .  (Of  course,  these  smoothness  conditions  can  be  weakened.) 


(6.9) 


(6.10) 


particular,  Q.  JRm  is  bounded, 
ators  with  respect  to  x  . ) 

We  prove  estimates  of  the 

v(x,t)  €  4i  (x,t)G 

and 

v(x,t)  €  *£(x,t)G£ 


In  Section  8.2  estimates  of  the  form 

v ( t )  €  4<  (t)  G  for  0  <  t  <  T  (6.11) 

are  treated,  where 

v(t)  €  Co (n)  is  defined  by  (v(t))(x)  =  v(x,t)  , 

Co (n)  =  (C0(n))n  is  the  set  of  continuous  functions  on  fi 

with  values  in  2Rn  , 

*  €  CQ  [0  ,T  ]  n  C.,(0,T]  , 

G  c  Co(ft)  satisfies  conditions  similar  to  those  described 
above  for  G  c  Hn  . 


In  the  basic  results  for  estimates  (6.11)  of  v  £  R  the 
boundary  terms  (6.7)  do  not  occur.  However,  in  verifying  the 
assumptions,  in  general,  certain  boundary  conditions  are  needed  for 
all  components  vi  for  which  the  coefficient  d^  does  not  vanish. 
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In  particular,  we  discuss  estimates  of  the  form 

(/  vT(x,t)  H  (x)  v (x , t )  dx)^  <  4.  (t)  (0  <  t  <  T)  (6.12) 

a 

where  H  (x)  =  (hjjx)  6^j)  with  h^  €  C0(fi)  ,  Imfx)  >0  on  f!  . 

This  estimate  is  equivalent  to  (6.11)  with 

G  =  {y  %  C"(fl)  :  W(y)  <  1} 

,  .  T  (6'13) 

w(y)  =  <(y,y )>  ,  <y ,n)>  =  J  y  (x)  H(x)  n(x)  dx  . 

n 

In  particular,  one  may  choose  H  (x)  =  I  or  H(x)  =  (d  1  (x)  6^)  , 

assuming  in  the  latter  case  that  all  coefficients  d.  are  strictly 
positive. 

The  results  in  [4]  on  estimates  (6.11)  for  parabolic  operators 
are  derived  from  a  theory  on  corresponding  estimates  for  a  first  order 
differential  operator  in  a  Banach  space: 


r v’ (t)  -  F(t,v(t) ) 

ft?v(  t)  = 

for  O  <  t  <  T 

v  (0)  -  r 

for  t  =  0 

These  abstract  results  will  not  be  described  here. 

Analogously,  estimates  on  abstract  differential  operators 
of  the  second  order  can  sometimes  be  used  to  derive  estimates  for 
elliptic  differential  operators  (see  [5]). 


The  various  estimates  described  here  require  different 
assumptions  on  the  form  of  the  matrices  D(x)  and  a (x)  and  on 
the  coupling  of  the  function  f  ,  where  both  types  of  assumptions 
are  related  to  *'ach  other  (see  the  discussion  below)  . 


In  [4]  we  proved  also  estimates  which  are  obtained  when  the 
sets  G  and  G^  in  (6.9),  (6.10)  and  (6.11)  are  replaced  by  its 
interior.  The  corresponding  theory  is  similar,  but  somewhat  simpler. 


As  mentioned  above,  the  estimates  proved  yield  invariance 
statements,  by  choosing  constant  functions  <|>  (4>&)  .  Invariance 

statements  (of  type  (U)  and  (E) )  have  been  proved  for  elliptic 
operators  in  [6-12],  for  parabolic  operators  in  [7,8,9,13-16],  for 
abstract  differential  operators  in  [17-20].  In  [21]  certain 
estimates  are  proved  which  can  be  considered  as  special  cases  of 
(6.11)  with  (|>(t)  *  const.  .  The  paper  [22]  treats  estimates  (6.12) 
with  H  (x)  »  I  and  4>(t)  =  exp(-xt).  In  [23]  estimates  v(t)  €  Gt 
are  derived.  For  further  references  see  [2-4]. 


7.1  Shape-invariant  bounds  j 

For  deriving  estimates  (6.3)  we  have  to  assume  that  all  functions 
d^  are  equal  and  all  functions  ai  are  equal: 

d^  —  d  ,  a ^  -  ct  (i  =  1 , 2, . . . ,n)  .  (7.1) 

We  shall  explain  below,  that  this  inconvenient  assumption  is  not 
necessary  for  more  general  estimates  (6.4). 

First,  we  consider  the  following 


Special  case:  G  is  convex,  f  does  not  depend  on  ux  .  (7.2) 

Here  the  theory  leads  us  to  differential  inequalities  for  4.  ,  which 
involve  a  parameter  n  e  lRn  .  These  inequalities  have  the  following 
form  (recall  (6.5)): 

-  o)(  n )  d  A4  +  W'(n)f(x,(J»n)  ^  W'(n)Mv  (7.3) 

on  the  set  of  all  (x,n)  with  x  €  fi  ,  W(n)  =  1  . 

-<i)(n)a-~ j  +  W' (n)Bn  <|>  ^  w*(n)Mv  (7.4) 

on  the  set  of  all  (x,n)  with  x  6  3S2  ,  W(n)  =  1  . 


Inequalities  of  this  type  occur  in  results  of  both  types  (U)  I 
and  (E)  ,  which  we  have  proved  for  this  special  case.  1 

The  results  of  type  (U)  are  derived  by  means  of  the  continuity  ! 
principle  (Section  4)  using  a  family  of  functions  4<^  =  4<  +  z^ 

(0  <  X  <  «>)  .  One  obtains  differential  inequalities  for  all  these 

functions  4>.  .  In  general,  one  can  choose  functions  z.  of  the  ! 

A  A  j 

special  form  z^  =  Xz  with  a  strictly  positive  function  z  .  Also, 
one  may  split  each  differential  inequality  for  4^  into  a  condition 
for  4-  and  a  second  condition  involving  z^  .  Then,  for  z^  =  Xz  j 

one  obtains  the  differential  inequalities  for  <|>  described  above  j 

and,  in  addition,  the  following  inequalities:  j 


-  u)(n)  d  Az  +  W' (n)  [f  (x,  (4»+Xz)n )  -  f(x,4<n)]  >  0 
on  the  set  of  all  (X,x,n)  such  that  X  >  0  ,  x  €  fi  and 
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w(n)  =  1 


v  ( x )  =  U(x)  +  Xz(x)n  ; 


(7.5) 


-  to  C  n )  a  +  W'(n)Bn  z  >  O 

on  the  set  of  all  (x,n)  with  x  £  9ft  ,  W(n)  =  1  . 

The  idea  is  to  prove  the  existence  of  suitable  functions  z 
for  a  certain  class  of  operators  M  .  Then,  for  each  operator  M 
in  this  class  the  following  statement  holds: 


If  the  inequalities  (7.3),  (7.4)  are  satisfied,  then 
v(x)  €  4>(x)G  for  x  £  fi  .  In  particular,  each  solution  v  of 
Mu  =  0  satisfies  this  estimate,  if  (7.3),  (7.4)  hold  with  Mv 
replaced  by  0  . 


In  [3]  the  application  of  these  results  is  explained  in  detail, 
in  particular,  the  use  of  side  conditions  on  v  such  as  the  one 
in  (7.5). 

Results  of  type  (E)  have  been  proved  for  the  special  case 
(7.2)  under  additional  assumptions  on  d  ,  o  ,  B  ,  f  ,  an  ,  such  as 
d(x)  >  d0  >  O  and  smoothness  conditions.  For  this  case,  the 
statements  have  the  following  form: 

If  the  inequalities  (7.3),  (7.4)  are  satisfied  with  Mv 

replaced  by  O  ,  then  the  boundary  value  problem  Mu  =  O  has  a 
solution  v  such  that  v(x)  €  4>(x)G  for  x  £  fi  . 


For  the  general  case,  where  f  may  depend  on  ux  and  G 
need  not  be  convex,  we  have  proved  results  of  type  (U)  .  Here,  the 
differential  inequalities  on  n  become  more  complicated.  In 
particular,  these  inequalities  involve  two  parameters  n  £  lRn  and 
q  £  ]Rn»m  .  For  example,  the  condition  for  4,  which  generalizes  (7.3) 
has  the  following  form,  where  Q(x,n,q)  is  the  trace  of  the  m*m 
matrix  d(x)qTw"(n)q  : 

-  w( n )  d  i(|i  +  Q (x, n #q) <l>  +  W' (n)f  (x,4,nf*q+n<|»x)  £  W'  (n)Mv  (7.6' 

on  the  set  of  all  (x,n,q)  with 

x  £  n  ,  w(n)  =  1  ,  w' (n)q  =  0  . 

Observe  that  the  parameter  q  is  not  bounded,  so  that 
assumptions  of  this  type  may  impose  strong  conditions  on  the  depend¬ 
ence  of  f(x,u,ux)  on  ux  .  We  shall  make  some  brief  remarks  concern- 


ing  conditions  on  the  coupling  of  f  with  respect  to  ux  . 

If  the  set  G  is  "strictly  convex"  and  d (x)  >  d0  >  0  ,  then 

the  term  Q(x,n,'g)  in  (7.6)  is  strictly  positive  (for  W' (n)q  =  0  ) 

and  W(y)  depends  on  all  components  yi  (i=1 ,2, . . . ,n)  .  Thus,  (7.6) 
constitutes  a  certain  quadratic  growth  condition  on  f  as  a  function 
of  ux  ,  but  f  may  be  strongly  coupled. 

If  W  depends  only  on  some  of  the  variables  y^  ,  say  y^ 

with  i  6  J  ,  then  Q(x,n,q)  does  not  depend  on  elements  q^  of 

q  with  j  €  J  .  Consequently,  the  condition  (7.6)  cannot  be  satis¬ 
fied  for  all  q  (with  W' (r  '  •  -  0  ) ,  if  f  depends  on  any  derivative 
SUj/Dx^  with  j  t  I  ,  \  for  very  special  nonlinear  functions 

f  . 

In  particular,  we  investigated  linear  operators  M  with 

Dirichlet  boundary  t Cixos  {or  the  case  W  (yj  =  yTy  .  These  statements 
on  linear  operators  contain  results  on  boundary  maximum  principles, 
more  general  than  those  of  [24,25,26], 


7.2  Estimates  v  £  fl  4^G^ 

For  estimates  of  the  more  general  form  (6.4),  where  Gjj.  is 
described  by  a  function  W£  ,  condition  (7.1)  can  be  replaced  by 
the  following  weaker 

Assumption  A  :  If  there  exists  an  Z  such  that  W£ (y) 

depends  on  yi  and  yj  ,  then 

d^  =  dj  and 

Example  1 ;  If  all  d^  are  pairwise  different,  then  each  W£ 
must  depend  on  one  variable  yj  only.  This  is  the  case,  for 
instance,  if 

N  =  2n  ,  Wi(y)  =  yi  ,  wn+i(y)  =  ~  Y±  (i=1 »2, . . . ,n) 

The  corresponding  statement  (6.4)  is  an  estimate  by  two-sided  bounds. 

Example  2:  In  order  to  obtain  an  estimate 

(v*  (x)  +  v^(x))^  <  ,  ( x )  ,  ( V3  (x)  +  vjj(x))2'  <  41 2  ( x ) 

in  the  case  n  =  4  ,  one  may  choose 

W,  (y)  =  y,  +  y“  ,  w2(y)  =  y^  +  y*  . 

Here,  we  need  to  require  that  d.,  =  d2  ,  d3  =  d4 
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Thus,  the  choice  of  functions  W£  which  do  not  depend  on  all 
yi  ,  has  certain  advantages;  but  it  has  disadvantages,  too.  Such  a 
choice,  in  general,  requires  stronger  conditions  on  the  coupling  of 
f(x,u,ux)  with  respect  to  ux 

In  Example  1 ,  each  component  f^  of  f  must  not  depend  on  any 
derivative  of  any  function  Uj  with  j  *  i  ,  except  for  very  special 
nonlinear  functions  f^  .  In  Example  2,  we  must,  in  general,  assume 
that  fi  and  f2  do  not  depend  on  any  derivative  of  u3  and  u«  , 
and  f3  and  f<*  do  not  depend  on  any  derivative  of  ui  and  u2  . 

For  a  better  understanding  of  the  conditions  required  see  the 
remarks  concerning  the  coupling  of  f  which  follow  formula  (7.6). 

For  the  more  general  estimates  considered  here  one  obtains  differen¬ 
tial  inequalities  for  the  functions  <j»£  analogous  to  (7.6),  with 
W  replaced  by  W£  ,  etc.  . 

We  proved  in  [3)  results  of  type  (U)  and  (E)  for  the 
special  case  (7.2),  and  results  of  type  (U)  for  the  general  case. 

As  an  example ,  we  derived  in  [3]  existence  statements  and 
numerical  error  estimates  for  the  steady  state  of  a  simple  reaction- 
diffusion  problem,  where 

Kx:  =  n£  =  i  ^£^G£  is  a  non~rectanc1ular  quadrangle  . 

This  choice  of  the  sets  G £  is  natural  because  of  the  special  struc¬ 
ture  of  such  problems. 


8.  Parabolic  operators 
8.1  Shape-invariant  bounds 

The  theory  in  Section  7  on  results  of  type  (U)  can  also  be 
applied  to  parabolic  operators,  as  discussed  in  [3].  One  obtains 
somewhat  stronger  (and  additional)  statements,  however,  by  the 
approach  described  in  [4],  We  shall  present  some  of  the  pertinent 
results  in  more  detail,  since  they  have  not  been  published  in  this 
form.  V.’e  use  the  definitions  and  notation  in  Section  6,  and  assume 
in  this  section  that  (7.1)  holds. 


Theorem  A  .  The  estimate 

v(x,t)  €  <|>(x,t)G  for  x€fi,  0  <  t  <  T  (8.1) 

holds  under  the  following  assumptions  (I)  ,  (II)  ; 
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(i)  u(n)(4't"d  A<i> )  +  Q  (x ,  n  , q)  +  W’  (n)  f  (x,  t,4>n ,4q+n4x) 

>  W‘ (n)Mv 


on  the  set 

of  all 

(x,t,n,q)  with 

x  e  n 

,  0  <  t  <  T  , 

W(n)  =1  ,  W*  ( n ) q  =  0  ; 

(8.2) 

(ii) 

o>(  n  H 

+  W'  (n )  Bn  4,  > 

W’  (n)Mv 

on  the  set 

of  all 

(x,t,n)  with 

x  £  3f2 

o 

A 

n- 

IA 

•«* 

W(n)  =  1  ; 

(8.3) 

(iii) 

v  (x,0) 

£  4<(x,0)G  for 

x  £  n 

(II)  Suppose  that  there  exists  a  function  z  £  S  such  that 
z(x,t)  >  O  on  n  *  [0,Tj  and 

(i)  u(n)(z  -  d  Az)  +  Q(x,n,q)z 

+  y  W' (n) [f (x,t/v,vx)  -  f (x,t,v-Xzn,vx-Azq-Xnzx) ]  >0 
on  the  set  of  all  (x,t,A,n»q)  which  satisfy  (8.2),  O  <  A  <  e  and 
v(x)  -  Az(x)n  =  <Hx)  ,  v ( x)  -  A z  (x) q  -  Xnzx(x)  =  iPx(x) 

(8.4) 

(ii)  ai(n)zv  +  W' (n)Bn  z  >  0 

on  the  set  of  all  (x,t,n)  which  satisfy  (8.3). 

Assumption  (II)  can  be  verified  under  rather  mild  conditions. 
For  example,  if  G  is  bounded,  the  following  statements  (a)  ,  (bj 
hold,  where 

q,p  €  ]Rn'm  and  ||p||  =  (Iik  (Pik)2)^  : 

(a)  Assumption  (II)  (i)  is  satisfied,  if  there  exist 
constants  yi  >  0  ,  c0  >  0  ,  c-,  ,  ca  ,  such  that 

u>(n)  s  pi  / 

Q(x,n,q)  >  c0||q||2  , 

W*  (n)  [f  (x,v,vx)  -  f  (x,v-«n,vx-p)  ]  >  c,a  -  ca(|p||  (8.5) 

on  the  set  of  all  (x, t , n ,q , a ,p)  which  satisfy  (8.2)»  a  £  0  and 
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v(x)  -  an  =  4-{x)n  ,  vx(x)  ~  P  =  td^tx)  ,  (8.6) 

and  if,  in  addition,  there  exist  a  function  <p  €  C2(n)  and  a 
constant  x  such  that 

<p  >  0  on  n,-dAq>  +  xip>0  on  ft  . 


(b)  Assumption  (II)  (ii)  is  satisfied,  if 
a  -|J  +  w'  (n)Bn  <p  >  0 

on  the  set  of  all  (x,ri)  with  x  £  3!)  ,  W(n)  =  1  • 


(One  verifies  by  formal  calculations,  that  under  the  assumptions 
made  above,  the  function  z(x,t)  =  <p  (x)  exp(Nt)  satisfies  (II)  (i) 
and  (ii)  ,  respectively,  if  N  is  sufficiently  large.) 


Remarks .  1)  Assumption  (8.5)  essentially  is  a  local  Lipschitz 
condition  on  f  ,  due  to  the  side  condition  (8.6).  Observe  also  that 
the  constants  Ci  and  c2  ,  which  may  depend  on  v  ,  need  not  be 
known  numerically. 

2)  In  the  case  of  a  Dirichlet  boundary  operator 

Mv(x,t)  =  v(x,t)  for  x£3il,0<t5T  , 

one  may  choose  <p  (x)  =  1  .  In  may  other  cases  a  function  of  the  form 
<p  (x)  =  rg  -  r2  can  be  used,  where  r  is  the  distance  of  x  from 
some  point  in  ft  and  r0  is  sufficiently  large. 

3)  If  Assumption  (II)  is  satisfied,  then  the  estimate  (8.1) 
follows  from  the  conditions  on  4>  described  in  Assumption  (I)  . 
Assumption  (I) (i)  constitutes  conditions  on  the  dependence  of 
f(x,t,u,u  )  on  u  ,  as  described  for  elliptic  operators. 

4)  In  [4]  we  explained  in  more  detail,  in  which  way  the  results 
can  be  used  for  an  error  estimation,  if  f(x,t,u,ux)  satisfies  a 
certain  quadratic  growth  condition  with  respect  to  ux  .  (See 
Corollary  l.c  in  [4]  and  the  text  succeeding  this  corollary.)  In 
particular,  one  can  obtain  in  this  way  estimates  (8.1)  with 

<|>(x,t)  =  «(x)  +  <p(x)  exp  (-pt)  , 

which  yield  statements  on  the  stability  of  solutions. 


In  some  of  the  results  in  [4],  we  did  not  use  the  side  condi¬ 
tions  (8.6).  For  this  reason,  we  formulated  there  an  assumption  which 
is  obtained  from  (8.5),  when  Cia  -  c2||p||  is  replaced  by 
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Ci a  -  c 2! |p ||  -  c 3||p ||  2  .  (See  Corollary  1b.) 


8.2  More  general  estimates 


In  [4]  we  described  also  results  of  type  (U)  on  estimates 


v(x,t)  £  <J,^(x,t)G£  (£=1 ,2, . . . ,N) 


(8.7) 


for  parabolic  operators.  For  such  estimates  condition  (7.1)  can  be 
replaced  by  Assumption  (A)  in  Section  7.2. 


In  addition,  [4]  contains  results  of  type  (U)  on  estimates 
of  the  form  (6.11): 


v ( t )  e  (Kt)G  (o  <  t  <  T)  , 


(8.8) 


where  (6.12)  is  a  special  case.  These  estimates  have  the  following 
important  advantage  to  estimates  of  the  form  (8.1)  or  (8.7): 

The  coefficients  d^  (and  a±  )  may  all  be  different  and,  at 
the  same  time,  there  is  no  restriction  on  the  coupling  of  f  . 


The  basic  results  on  estimates  (8.8)  contain  conditions  on 
certain  function  z  (or  a  family  of  functions  z\  )  which  are 
analogous  to  those  in  Theorem  A  .  We  shall  not  describe  this 
general  theory  here.  Instead,  we  formulate  a  more  special  result 
on  estimates  (6.12). 

Here,  we  define  ^  ,  )as  in  (6.13)  and  use  abbreviations  such 
as 

</),  f  (v,vx))  (t)  =  J  nT  (x)  f  (x,t,v(x,t)  ,vx  (x  ,t) )  dx 

Moreover,  for  sufficiently  smooth  n  £  Cq(H)  let 

l(n)  =  /  nT(x)  H  (x)  (-  D  (x)  An(x))  dx  (8.9) 

n 

and  denote  by  ||  njl  some  semi-norm  of  the  matrix  nx  (see  the 
example  below) . 


Theorem  B.  The  estimate 

(/  vT(x,t)  H(x)  v(x,t)  dx)^"  <  th  (t)  (0  £  t  <  T) 

fi 

holds  under  the  following  assumptions  (i) ,  (ii)  ,  (iii): 

(i)  Suppose  there  exist  constants  6  >  0  ,  vo  »  Pi  »  Pa  >  0 
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Ci  ,  c2  >  O  such  that 

Kn)  ^  Usll^l!2  +  PillnJI  +  yo  (8.10) 

and 

<^ri,f(v,vx)  -  f  (v-an  ,  vx~arlx)/>  (t)  >  C-i  a  -  c2o||r^||  (8.11) 

for  all 

a  e  (0,6]  ,  t  e  (o , t ]  ,  n  6  c"(n)  (8.12) 

which  satisfy 

^n,n^  =  1  and  v(x,t)  =  ( 4,  ( t )  +  a  )n(x)  for  x  €  n  .  (8.13) 

(ii)  Suppose  that 

<l>t(t)  +  l(n)4»(t)  +  ^n,  f  (<]>n ,  4>nx)^  (t)  >  (n,Mv)  (t)  (8.14) 

under  the  side  conditions  (8.12),  (8.13). 

(iii)  <(v,v^>  (0)  <  4, 2  (O) 

Remarks .  1)  Due  to  (8.13)  the  function  n  occuring  in  (1) 
and  (ii)  is  sufficiently  smooth. 

2)  Estimates  of  the  form  (8.10),  in  general,  can  be 
obtained,  if  v  satisfies  certain  boundary  condition.  One  uses 
partial  integration  and  the  eigenvalue  theory  of  elliptic  differential 
operators. 

3)  The  assumption  (8.11),  (8.13)  describes  a  certain  local 

Lipschitz  condition,  which  is  satisfied  for  sufficiently  smooth  f  . 

4)  In  the  differential  inequalities  (8.14)  for  4  ,  one 
will  use  an  estimate  of  the  form  (8.10)  where,  however,  p2  need 
not  satisfy  p2  >  0  .  If  p2  >  0  ,  then  (8.14)  constitutes  a  certain 
quadratic  growth  condition  on  f  as  a  function  of  ux  . 

Example :  Suppose  that 

n  =  (x  €  ]Rn  :  |xi|  <  1  for  all  i  }  ,  H  (x)  =  D(x)  =  I  , 

and 

v(x,t)  =  0  for  x  E  9S1 

Then 

Kn)  =l!^ll2  with  ||  ||  =  (J  (x)  dx)^ 
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and 


on  u 
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III.  NUMERICAL  APPLICATIONS 


The  results  described  above  have  many  applications  in  numerical 
analysis.  For  example,  they  can  be  used  as  tools  to  prove  the 
convergence  of  iterative  procedures  or  the  convergence  of  difference 
methods  for  differential  equations.  This  was  described  for  linear 
problems  in  [1]  (see  Section  2  above).  The  application  to  error 
estimations  is  discussed  in  the  following  section. 


9.  Approximation  methods  and  a  posteriori  error  estimates  for 
two-point  boundary  value  problems 


9.1  Error  estimates 


In  principle,  each  of  the  estimation  methods  described  in  the 
previous  sections  can  be  applied  to  obtain  a  posteriori  error  bounds 
for  approximate  solutions  of  equations.  (For  the  application  of 
inverse-monotonicity  see  [27,28].)  The  question  arises,  whether 
one  can  construct  algorithms  of  error  estimation  which  on  one  hand 
yield  sufficiently  good  error  bounds  and,  on  the  other  hand,  do  not 
require  too  much  additional  work.  We  investigated  these  possibilities 
for  the  case  of  (scalar-valued  and  vector-valued)  two-point  boundary 
value  problems  of  the  second  order,  as  part  of  the  research  project 
described  in  this  report.  (For  methods  of  error  estimation  for 
initial  value  problems  see  [29,30].)  The  work  on  boundary  value 
problems  has  not  yet  been  completed.  Nevertheless,  it  seems 
worthwhile  to  report  briefly  on  our  experiences. 


Let  us  first  describe  the  general  principle  of  error  estimation. 
For  a  given  equation 
Mu  =  r 


denote  by 


The  error 


u* 

a  solution. 

w 

an  approximation  solution. 

V  =  u*  -  w 

the  error  of  w 

v  satisfies 

the  equation 

Mv  =  d 


with  the  known 


27 


defect  d  =  d[w]  =  -  Mw  +  r  ,  and  Mv  =  M(w+v)  -  Mw 

After  this  reformulation  of  the  given  problem  one  tries  to  apply 
one  of  the  estimation  methods  described  in' the  previous  sections, 
in  order  to  obtain  an  inclusion  statement  v  €  K  for  the  error. 


Some  first  results  on  the  error  estimation  for  two-point  bound¬ 
ary  value  problems  are  described  in  [1,  Sections  11.5.2  ,  III. 6. 3, 
and  Example  V.4.12].  The  example  mentioned  is  a  vector-valued 
problem  with  a  singular  differential  operator,  and  the  estimate  was 
one  by  pointwise  norm  bounds.  Here,  we  shall  restrict  ourselves  to 
considering  scalar-valued  equations  of  the  form 

-  u"  (x)  +  f(x,u(x))  =  0  for  0  <  x  <  1 

with  linear  Sturm-Liouville  boundary  conditions,  and  estimates 

|  v  (x)  |  =  |u*(x)  -  W  ( X )  j  <  (Mx)  (0  <  x  <  1) 

with  (Ji  €  C2[0,1]  .  Suppose,  for  simplicity,  that  the  approximate 
solution  w  £  C2  [0,1]  satisfies  the  boundary  conditions,  and  write 

d (x)  =  d[w] (x)  =  w" (x)  -  f(x,w(x)) 

For  continuous  f  (on  [0,1]  *  ]R  )  one  obtains  the  following 
statement: 


If  >  0  satisfies  the  homogeneous  boundary  conditions  and  if 

-  <{-"  +  f(x,w+4/)  -  f(x,w)  >  d  on  (0,1) 

and 

-  4."  +  f(x,w)  -  f(x,w-<|))  >  -  d  on  (0,1)  , 

then  there  exists  a  solution  u*  6  Ca[0,1]  such  that  |  u*  -  w|  <  4* 
on  [0, 1 ]  I 

The  application  of  this  result  poses  essentially  three  problems 

1)  The  computation  of  an  approximate  solution  w  . 

2)  The  estimation  of  the  defect  d[w]  on  the  entire 
interval  [0,1]  . 

3)  The  construction  of  a  function  4>  ,  which  solves  the 
differential  inequalities. 


To  1:  Approximation  methods  will  be  discussed  in  the 
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succeeding  section. 

To  2:  The  defect  d(x)  usually  is  a  function  of  a  very 
complicated  form.  One  has  to  replace  it  by  a  simpler  bound,  say  a 
constant  r  : 


|d(x) I  <  r  (0  <  X  <  1)  .  (9.1) 

In  the  general  case,  the  calculation  of  a  good  bound  r  ,  such  that 
this  inequality  holds  for  all  x  €  [0,1]  ,  is  the  most  difficult 
and  time-consuming  part  of  the  procedure. 

There  are  cases,  however,  where  such  an  estimate  is  comparibly 
easy.  If  d(x)  is  a  polynominal,  then  r  can  be  obtained  by 
calculating  finitely  many  values  d(xk)  .  Here,  one  uses  a  result  of 
[ 31 ] .  In  such  cases,  the  estimate  (9.1)  often  requires  much  less  work 
than  the  computation  of  w  . 

This  method  of  estimating  the  defect  can  also  be  used,  if  [0,1] 
is  divided  into  a  few  intervals  and  d(x)  is  a  polynomial  on  each 
of  these  subintervals. 

More  general  cases  are  still  being  investigated. 


To  3:  In  general,  this  is  the  easiest  part  of  the  procedure 
(see  [1,  Section  III. 6.1]). 


9.2  Approximation  methods 


For  carrying  out  the  error  estimation  we  need  to  calculate  an 
approximate  solution  w  ,  which  is  defined  for  all  x  €  [0,1]  . 

There  exist  several  methods  and  programs  to  obtain  such  approximations 
(see  [32,33],  for  example). 

In  order  to  investigate  the  power  of  the  estimation  methods, 
we  first  used  polynomial  approximations.  A  series  of  approximation 
methods  for  obtaining  such  (global)  polynomials  were  tested.  We 
finally  arrived  at  the  following 

Method  P  :  (Collocation  using  integrated  Legendre  polynomials 
and  Chebychev  collocation  points.)  Write  the  polynomial  w(x)  in  the 
form 


n 

w (x)  =  <p0  (x)  +  jTj  ai  <Pi(x)  , 

where  ?o(x)  is  linear  and  the  functions  <Pi(x)  are  once  (or  twice) 
integrated  Legendre  polynomials  (depending  on  the  boundary  conditions) « 
Calculate  the  constans  a^  such  that 


_  J 
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d[w](xk)  =0  (k=1 ,2, . . . ,m)  , 

where  denote  the  zeros  of  the  m-th  Chebychev  polynomial 

(transformed  onto  [0,1]  ). 


We  obtained  excellent  results  (concerning  accuracy  and  computing 
time)  with  method  P  for  a  series  of  problems  whose  solutions  did  not 
behave  "too  badly".  For  example,  the  problem 

-  u"  -  [90  ctg  ■g-Trd+x)  +  60  tan  ^•■n(1+x)]u'  -  630  u  =  0  , 

u  (0)  =  0  ,  u(1)  =  5 

could  be  solved  without  any  difficulty.  The  solution  u*  increases 
very  fast  from  u*  (0)  =0  to  its  maximal  value  224  ,  which  is 
attained  at  approximately  x  =  0.022  . 

For  a  series  of  problems  of  this  type  we  compared  this  polynomial 
approximation  with  methods  of  spline  approximation,  using  programs 
of  de  Boor  [32]  and  Ascher  et  al.  [33].  It  turned  out  that  for  such 
problems  of  "medium  difficulty"  method  P  gave  by  far  the  best  results 
(see  [34]). 

Of  course,  this  does  not  mean,  that  one  should  use  polynomials 
instead  of  spline  functions  in  general  differential  equation  solvers. 
For  example,  method  P  failed  for  certain  singular  perturbation 
problems  with  very  small  perturbation  parameters,  while  the  other 
methods  still  were  applicable.  Our  results  suggest,  however,  to 
construct  programs  for  the  approximation  by  piecewise  polynomial 
functions  where  the  polynomials  in  the  subintervals  can  have  dif¬ 
ferent  orders  and  very  high  orders.  Then  one  may  not  need  many 
subintervals,  after  one  has  found  out,  where  the  boundary  layers, 
transitions  layers  etc.  occur.  Such  a  procedure  would  also  simplify 
the  a  posteriori  error  estimation. 
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